Collective modes of the massless Dirac plasma 
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We develop a theory for the long-wavelength plasma oscillation of a collection of charged massless 
Dirac particles in a solid, as occurring for example in doped graphene layers, interacting via the 
long-range Coulomb interaction. We find that the long-wavelength plasmon frequency in such a 
doped massless Dirac plasma is explicitly non-classical in all dimensions with the plasma frequency 
being proportional to We also show that the long wavelength plasma frequency of the D- 

dimensional superlattice made from such a plasma does not agree with the corresponding D -I- 1 
dimensional bulk plasmon frequency. We compare and contrast such Dirac plasmons with the well- 
studied regular palsmons in metals and doped semiconductors which manifest the usual classical 
long wavelength plasma oscillation. 
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A collection of charged particles (i.e. a plasma), elec- 
trons or holes or ions, is characterized by a collective 
mode associated with the self-sustaining in-phase den- 
sity oscillations of all the particles due to the restoring 
force arising from the long-range 1/r Coulomb poten- 
tial. The classical plasma frequency in three-dimensional 
(3D) plasmas [Ij] is well-known to be W3 = {ATrnse^ /rriY^^ 
where e and m are respectively the charge and the mass 
of each particle, and ns is the 3D particle density. (In 
this paper, we use uju and tid as the Z)-dimensional 
long-wavelength plasma frequency and particle density 
respectively.) A solid state degenerate plasma 0, 0] 
exists in metals and doped semiconductors where free 
carriers can move around quantum mechanically in the 
ionic lattice background. Such a degenerate quantum 
plasma has the quantized version of exactly the same 
collective mode, the so-called plasmon 0, S I3|' which 
dominates the spectral weight of the long-wavelength el- 
ementary excitation spectrum of an electron liquid. (We 
will use the world 'electron' generally throughout this 
paper to indicate either electron or hole.) The collec- 
tive plasmon modes of solid state quantum plasmas have 
been extensively studied experimentally and theoretically 
over the last sixty years in both metals and doped semi- 
conductors. In the present work, we study theoretically 
the collective plasmon mode in a solid state plasma of 
massless Dirac fermions, as occurring for example, in 2D 
graphene layers. We define the Dirac plasma as a system 
of charged carriers whose energy-momentum dispersion 
is linear, obeying the Dirac equation. 

Our main qualitative result is that the massless Dirac 
plasma is manifestly quantum, and does not have a classi- 
cal limit in the form of an ?i-independent long-wavelength 
plasma frequency, in a striking contrast to the corre- 
sponding parabolic dispersion electron liquids familiar 
from the extensive study of plasmons in metals and semi- 
conductors 0, 0, HI. The long- wavelength plasmon fre- 
quency of a Dirac plasma is necessarily quantum with 



^fV appearing manifestly in the long wavelength plasma 
frequency in D = 1, 2, 3 dimension (and in between). 
By contrast the long wavelength plasma frequency of or- 
dinary electron liquids is classical, and quantum effects 
show up only as nonlocal corrections in higher order wave 
vector dispersion of the plasmon mode. This is quite 
unexpected in view of the popular belief that the long 
wavelength quantum plasmon dispersion is necessarily a 
classical plasma frequency 0, 13, [j] . The popular belief 
seems to be true for the usual parabolic energy disper- 
sion, but not for the linear Dirac spectrum. 

We start from the fundamental many-body formula 
defining the collective plasmon mode in an electron sys- 
tem: 



e{q,Lo) = l~viq)niq,co)=0, 



(1) 



where e{q,uj) is the wave vector (q) and frequency (lo) 
dependent dynamical dielectric function of the system, 
with Il{q,uj) the irreducible polarizability and v{q) the 
Coulomb interaction between the electrons in the wave 
vector space. The zero of the dielectric function in Eq. (1) 
signifies a self-sustaining collective mode, with the solu- 
tion of Eq. (1) giving the plasmon frequency as a function 
of wave vector. We first recapitulate the known results 
for the parabolic dispersion electron system before dis- 
cussing the novel collective dispersion for massless Dirac 
plasma. 

The Coulomb interaction in the wave vector space is 
given by the appropriate _D-dimensional Fourier trans- 
form of the Coulomb interaction v{r) — e^/rer 
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where we have introduced a background dielectric con- 
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stant {k) which, in general, differs from unity in semi- 
conductor based electron systems, and Kq is the zeroth- 
order modified Bessel function of the second kind. We 
note that Kq{x) ~ | ln(x)| for x ^ 0, and the length 'a' 
in the ID Coulomb interaction in Eq. (2c) characterizes 
the typical lateral confinement size of the ID electron 
system (ES) which is obviously necessary in defining a 
IDES. 

The irreducible polarizability function Il{q,uj) of an 
interacting ES is, in general, unknown since self-energy 
and vertex corrections cannot be calculated exactly. A 
great simplification, however, occurs in the long wave- 
length limit {q — > 0) when the dielectric function, and 
consequently, the plasmon frequency is determined en- 
tirely by the noninteracting irreducible polarizability, the 
electron-hole 'bubble' diagram. The noninteracting irre- 
ducible polarizability is given by the expression: 



U{q,uj) 



{27t)D hoj + & - 6+9 



F{k,q), (3) 



where 6 is the single particle energy dispersion i.e. 
6 — fi^k'^/2m for parabolic systems (and 6 = hvpk 
for the massless Dirac plasma), is the Fermi distri- 
bution function, and F{k, q) is the overlap form factor 
due to chirality. For non-chiral systems F{q, k) = 1. The 
factor '(7' in Eq. (3) is the degeneracy factor: g — g^gy 
where g^ (=2) is the spin degeneracy and g^ is the valley 
or pseudospin degeneracy. 

Putting 6 = fi^k'^/2m, we can easily calculate Eq. (3) 
upto the leading order in wave vector (i.e. the long wave- 
length limit) to obtain 
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Combining Eqs. (l)-(4) we immediately obtain the 
well-known long-wavelength plasma frequency in a D- 
dimensional ES: 
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where denotes the long-wavelength {q —> 0) plas- 
mon mode in the D-dimensional parabolic dispersion ES 
(with the carrier density n d per unit D-dimensional vol- 
ume) where the one particle energy is given by ^ = 
fi^fc^/2m /2m ~ mv^ /2 classically (where p — hk \s 
momentum). The long wavelength plasmon frequencies 
for parabolic dispersion systems given in Eq. (5) are, of 
course, well-known and have been verified experimentally 
extensively 0, 0, Our purpose of deriving Eq. (5) is 
the explicit demonstration, to be contrasted below with 



the corresponding massless Dirac plasma, that the long 
wavelength plasmon frequency w^P-* for parabolic systems 
is completely classical since 'fi' does not appear in the 
leading term of Eq. (5) in any dimension. (Note that the 
right hand side of Eq. (5) has the explicit dimensionality 
of time inverse, i.e. a frequency, in each dimension since 
no has the dimension of (length)"^ in I?-dimension.) 
The second order dispersion correction term in Eq. (5), 
i.e. the 0(g^, g^/^, g^) term in I? = 3, 2, 1 respectively, 
is fully quantum mechanical (i.e. '/i' shows up explicitly 
in the non-local wave vector corrections), and is affected 
by interaction corrections (both self energy and vertex 
corrections to the irreducible polarizability). 

Now we consider plasmons in the Z3-dimensional mass- 
less Dirac plasma, where the single-particle energy dis- 
persion is linear, i.e. 6 — ^^Flk] — > vp classically, in 
D = 1, 2, 3. The long wavelength quantum plasmon 
dispersion is still defined by the set of formula given 
by Eqs. (l)-(3) with the explicit form of the nonin- 
teracting irreducible polarizability being calculated with 
6 = hvF\k\ in Eq. (3). 

The long wavelength (g — > 0) form for the noninteract- 
ing irreducible polarizability (Eq. (3)) can be calculated 
for linear energy dispersion relation (i.e. 6 = hvpk) in 
all dimensions, giving 
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where kp is the Fermi momentum of the system and T{x) 
is the Gamma function. (We note that the chirality factor 
F{k, q) in Eq. (3) does not influence the long wavelength 
limit.) 

Combining Eqs. (1), (2), and (6), we get the following 
for the long wavelength plasmon frequency, , 
1, 2, 3 Dirac plasma: 
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(7b) 
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where we have introduced the dimensionless fine struc- 
ture constant rs(= /{kUvf)) for notational simplicity. 

Comparing Eqs. (5) and (7) we see that a;'-^-' and w^'^ 
have one important similarity and several striking differ- 
ences. The similarity is that the plasmon dispersion (i.e. 
the power law dependence of the plasma frequency on 
wave vector) is the same in the parabolic system and the 
massless Dirac plasma for all D. This is indeed required 
under very general principles since for any Coulomb sys- 
tem, the long wavelength plasmon dispersion is set by 
the continuity equation (or equivalently, by particle con- 



servation) to be oJoiq 0) ^ g 



(3~-D)/2 



as is obeyed by 



both wj^^ and lo 
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The most striking qualitative feature of UjJ in Eq. (7), 
in sharp contrast with the usual in Eq. (5), is that 
appears explicitly in the leading term, not just the sub- 
leading nonlocal corrections. A simple dimensional anal- 
ysis of Eq. (7) shows that uj'f^ ~ 0{h~^/'^) in all dimen- 
sions in contrast to the 0{ff') purely classical behavior 
of Lj^^^ in Eq. (5). There is no classical plasma frequency 
in the massless Dirac plasma, i.e. the long-wavelength 
plasma frequency for ES with linear dispersion explicitly 
depends on and is therefore, by definition, nonclassi- 
cal. This absence of a classical long wavelength plasma 
frequency in the Dirac plasma is a direct manifestation 
of the relativistic Dirac nature of the underlying quan- 
tum description, and such a Dirac plasma does not have 
a classical plasma frequency. Note that the nonclassical 
nature of the long wavelength plasma oscillation of the 
Dirac plasma is independent of the chirality or gapless- 
ness of graphene, and arises primarily from the linear 
Dirac spectrum. 

Associated with the appearance of \/\/h in the long- 
wavelength plasma frequency of the Dirac plasma are sev- 
eral other interesting properties distinguishing it from 
the standard parabolic dispersion Schrodinger plasma: 
(1) The density dependence of the Dirac plasmon is dif- 
ferent from the regular plasmon 0, — iii particu- 
lar, the density dependence is weaker in the sense that 
Lo'"!^ cx v}/^, n^/^, in D — 3, 2, 1 respectively in 
contrast with c>c ^/n in all dimensions. In general, 
the plasmon frequency in the Dirac plasma is given by 
(jj^-* oc rS'^^^'^/'^^ . (2) The ID Dirac plasmon frequency 
is curiously density independent. (3) The quantum cou- 
pling parameter (i.e. the effective fine structure constant) 
shows up explicitly in the long- wavelength Dirac plasmon 
frequency, w^' cx ^Jf^. (4) The long wavelength Dirac 
plasmon goes as hr^^'^ for all dimensions whereas 
the long wavelength regular plasmon wj^' goes as m~^/^ 
in all dimensions. 

Before concluding, we consider another interesting and 
peculiar feature of the Dirac plasmon distinguishing it 
from the regular plasmon. We consider collective modes 
of periodic arrays of 2D Dirac plasma layers (for exam- 
ple, a graphene superlattice made of parallel 2D graphene 
sheets in the direction transverse to the 2D graphene 
plane) and of ID Dirac plasma nanoribbons (i.e. a 
graphene superlattice made of identical ID graphene 
nanoribbons placed parallel to each other in the 2D 
plane). Collective plasmon modes of such 2D [9|] and 
ID [lOj superlattices have been theoretically studied in 
the context of regular parabolic systems, and have been 
experimentally observed in doped GaAs multi-quantum 
well and multi-quantum wire structures. 

The main physics to be considered in describing the 
collective plasmon modes of such superlattices is the in- 
clusion of the inter-layer or inter-ribbon Coulomb inter- 



action, which will necessarily couple all the layers (or the 
ribbons) due to the long range nature of the Coulomb 
potential. This changes the fundamental collective mode 
equation (Eq (1)) to an infinite matrix equation: 

\5u' -Vu'{q,u^)I^v{q,u^)\=Q, (8) 

where 11; = 11 is the irreducible polarizability of each 2D 
layer (or ID ribbon), which is exactly the same polariz- 
ability considered in Eqs. (3) and (4). In Eq. (8), vw is 
the Coulomb interaction between the I and the /' layer 
or ribbon in the periodic array, which is given by 

vu. = ^e-«'^l'-''l D = 2, (9a) 
nq 

2e2 

vw = —[KQ{qa)+Ko{qd\l-l'\)] D = l, (9b) 

where 'd' is the superlattice period (to be distinguished 
from the length 'a' in D = 1 which defines the lateral 
width of each ribbon). 

The periodic invariance of the superlattice and the as- 
sociated Bloch's theorem allow an immediate solution of 
the infinite-dimensional determinantal equation defined 
by Eq. (8), leading to the following collective plasmon 
bands for the superlattice structure: 

Cj2s{q;k) = uj2{q)S2{q,k) D ^ 2, (10a) 
u;is{q;k) ^ u;i{q)Si{q,k) D = 1, (10b) 

where lojjs is the plasmon band frequency for the su- 
perlattice {D = 2 for the multilayer and = 1 for the 
multiribbon periodic arrays) and uju is the corresponding 
2D {D = 2) and ID {D = 1) plasmon modes discussed 
in Eqs. (5) and (7). The wave vector q in Eq. (10) is the 
same conserved 2D or ID wave vector in each individual 
2D layer or ID nanoribbon defining the plasmon disper- 
sion relation uJoiq) whereas the additional wave vector 
fc is a new continuous parameter defining the superlat- 
tice plasmon band (arising from the periodicity in the 
array structure). The band wave vector k is restricted 
to the first superlattice Brillouin zone, k < ir/d, in the 
reduced zone scheme. For the 2D layer superlattice, if 
each layer is assumed to lie in the x-y plane, then k = qz 
is along the superlattice direction of the z-axis. For the 
ID ribbon superlattice, if each ribbon is assumed to be 
along the x-axis (i.e. q = qx) with a width of 'a' defining 
the ribbon in the y-direction, then k = qy is along the 
superlattice direction of the y-axis. 

The function So in Eq. (10) is a form factor arising 
from the Coulomb coupling between all the layers and 
the ribbons forming the periodic array, and is given by 

52 = ^e-<zl^-^'l<i-^9.1^-^'M ^ ^ (11a) 
I' 

Si = [Ko{q\l - l'\d) cos{lqyd) + Ka{qa)] D = I. (lib) 
I' 
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Combining the above equations for superlattice plas- 
mons, we get the following long wavelength {q 0) plas- 
mon bands for 2D and ID arrays in parabolic and Dirac 
plasma systems, respectively: 

- ,,(P,^)^„^^ sinh(qd) "1 



2s 



^(q) 



.{p,i) 



cosh{qd) — C08{qzd) 



Kciiqa) + 2 KQ{nqd) cos{qynd) 



1/2 



(12a) 



(12b) 



Eq. (12) above defines plasmon bands for superlattice 
arrays made out of periodic 2D layers and ID ribbons in 
parabolic and linear plasma systems. 

An interesting quantum feature of uj^^^^{q,k) is appar- 
ent when one looks at the long-wavelength plasmon {q — > 
0) at the band edge fc = 0, and compares a)^g(g, fc = 0) 
with Q^j^l{q,k^ 0). We get 

o2\ 1/2 



^^i\q;q.^O) = 



Kin 
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Km 
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with = —, (13a) 
d 

2 

with n2 — (13b) 



and 



'^2! (9; Qz = 0) = 



1/4 




VF, (14a) 



(14b) 



We note that Eq. (13) for the usual parabolic electron 
plasma has the appropriate physical limit at the band- 
edge A: = 0, where the D-dimensional superlattice plas- 
mon should have the precise character of the correspond- 
ing [D + l)-dimensional bulk plasmon in the long wave- 
length limit, and indeed uj^J{q, k = 0) and oj^J{q, k — 0) 
are identical to the corresponding 3D and 2D plasmons 
(in Eq. (5)) respectively with — n^ja and ni = ni/a. 
This is exactly what one expects since the D-dimensional 
superlattice "loses" its discrete periodic structure for 
k = Q and simply becomes the {D + l)-dimensional reg- 
ular plasmon at long wavelength. 

Amazingly this does not happen for the Dirac plasma, 
where the _D-dimensional superlattice plasmon for k = 
does not become the corresponding {D ^ 1) -dimensional 
bulk plasma frequency as one expects intuitively. In par- 
ticular, CL!2g{q, k = 0) would agree with the corresponding 
3D Dirac plasmon CL]^^\q) (in Eq. (7)) only if we define 
the corresponding effective 3D density to be 

n, = {97rg/m)'/^{n2/d'f\ (15) 

rather than the intuitive definition = n2/d. For the 
ID superlattice Dirac plasmon, the situation is quali- 
tatively different since (q, k — 0) does not depend 



at all on the carrier density, and only the following 
strange substitution provides a correspondence between 
co^/^,iq,k = 0) and (9^0): 



fi2 = 9/{'^(P), 



(16) 



which is a constant for all carrier density. This strange 
lack of correspondence between Cj'"j^^{k — 0) and is 
again a manifestation of the peculiar quantum nature of 
the Dirac plasma, where the band-edge plasmon at fc = 
in the Z?-dimensional superlattice does not reduce to the 
corresponding bulk plasmon in [D + l)-dimension, as it 
does for the ordinary parabolic ES. 

In summary, we have found that the long-wavelength 
plasma frequency of a masslcss Dirac plasma with lin- 
ear carrier energy dispersion is non-classical with an ex- 
plicit 1 / \/h appearing in the plasma frequency. This is in 
sharp contrast with the widespread expectation that the 
long- wavelength plasmon is a classical plasma oscillation 
- in fact, a massless Dirac plasma has no classical anal- 
ogy. In addition, the long-wavelength Dirac plasma fre- 
quency depends explicitly on the coupling constant ( "the 
fine structure constant"). We have also shown that the 
long wavelength plasma mode of a i'-dimensional super- 
lattice of massless Dirac plasma does not reduce to the 
corresponding [D -f l)-dimcnsional bulk plasmon, as one 
would have expected intuitively. All of these peculiar re- 
sults follow from the fact that a massless Dirac plasma is 
fundamentally non-classical since the energy dispersion 
E — vp characterizing a system with constant velocity 
(but variable momentum) simply cannot happen in clas- 
sical physics. We believe that our predictions can be 
tested in doped graphene layers and multilayers, and in 
doped graphene ribbons and multiribbons arrays using 
electron scattering [11], light scattering [l3|, or infrared 
[ist spectroscopies. But the real importance of our re- 
sults is conceptual as we establish a strange quantum 
behavior in the graphene world of a Dirac plasma where 
the long wavelength plasmon is explicitly non-classical 
in contrast to plasmons in ordinary semiconductors and 
metals whose long wavelength limit is necessarily a clas- 
sical plasma frequency. 
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